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ABSTRACT 

Context. A small group of bipolar protostellar outflows display strong emission from shock-tracer molecules such as SiO and CH 3 OH, 
and are generally referred to as "chemically active." The best-studied outflow from this group is the one in L 1 157. 
Aims. We study the molecular emission from the bipolar outflow powered by the very young stellar object HH 114 MMS and compare 
its chemical composition with that of the L 1 157 outflow. 

Methods. We have used the IRAM 30m radio telescope to observe a number of transitions from CO, SiO, CH 3 OH, SO, CS, HCN, 
and HCO + toward the HH 114 MMS outflow. The observations consist of maps and a two-position molecular survey. 
Results. The HH 114 MMS outflow presents strong emission from a number of shock- tracer molecules that dominate the appearance 
of the maps around the central source. The abundance of these molecules is comparable to the abundance in L 1 157. 
Conclusions. The outflow from HH 1 14 MMS is a spectacular new case of a chemically active outflow. 



Key words. Stars: formation - ISM: abundances - 
lines: ISM 



Introduction 



ISM: jets and outflows - ISM: individual (HH 114 MMS) - ISM: molecules- Radio 



^ Bipolar outflows from young stellar objects ( YSOs) are routinely 
Identified by their mechanical impact on the surrounding cloud, 
^ which produces characteristic lobes of blue- and red-shifted ma- 
i_i terial. A small group of outflows also show evidence of chemical 
impact in the form of enhanced abundance of shock-sensitive 
,— ' species such as SiO and CH3OH. These outflows are often re- 
^ ferred to as "chemically active," and their study offers a unique 
opportunity to investigate both ambient g as acceleration and 
t ne co mplex reactions of shock chemistry (Tafalla & Bachiller 
^|20Tfl). 

_J ■ Here we report the serendipitous identification of a chem- 
^ ically active outflow driven by HH 114 MMS. HH 114 MMS 



is locat ed in the L 15 89 cloud, at an estimated distance of 
450 pc (Mathieu 2008). It was discovered at mm-wavelengths 
bv lChiniet all (119971). and was la ter mapped at cm wavelengths 
. ^ by Rodriguez & Reipurthl ( fl996h . Both HH 1 14 MMS and the 
k> neighboring source IRAS 05155+0707 lie near the center of the 
^ %6-pc long chain of H erbig-Haro (HH) objects HH 1 14/HH 1 15 
d dReipurth et al.l [l997l) . although the relation between the two 
YSOs and the HH chain remains unclear. Each YSO drives a 
separate CO outflow, and there is n o evidence for an y interac- 
tion between the two. As shown by iLee etal (2002), the blue 
lobe of the HH 1 14 MMS outflow approximately coincides with 
HH 1 14D, but the large-scale orientation of the HH 1 14/HH 1 15 
system differs significantly from the direction of the two CO out- 
flows. While it is possible that IRAS 05155+0707 drives the 
large-scale HH chain (if its initially E-W direction changes to 
SE-NW), it is unlikely that the HH 1 14 MMS outflow plays any 



* Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
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role in the large-scale HH 1 14/HH 1 15 system, since that would 

require a change of outflow direction close to 90 deg rees. 

HH 114 MMS is a highly embed ded object. IChini et alj 



dl997h suggested it is a class object (lAndre et al.l 



19931T and 

determined an upper limit to its luminosity of about 25 L Q . Fur- 
ther evidence for its extreme youth comes from the high collima- 
ti on of its molecular outfl ow shown by the interferometric maps 
of lArce & Sargent) (120061) . 



2. Observations 

The HH 1 14 MMS bipolar outflow was observed in the CO(2-l) 
line with the IRAM 30m telescope during December 2008 using 
the dual-polarization 3x3 HERA array operating in wobbler- 
switching mode with a throw of 240" and a frequency of 0.5 Hz. 
The signal from HERA was sent to the VESPA autocorrelator, 
which was configured to produce spectra with a resolution of 
1.25 MHz (1.6 km s 1 ) over a 280 MHz (365 km s 1 ) passband. 
These spectra were later processed with the CLASS softwar43and 
converted to the main beam brightness temperature scale. 

Additional IRAM 30m observations of the HH 1 14 MMS 
outflow were made in February 2011. An initial survey for se- 
lected molecular lines in the 1, 2, and 3 mm wavelength bands 
was carried out toward several positions. In view of the bright 
lines of shock-tracer species, additional positions were observed 
to map the full extent of the emission. These shock-tracer obser- 
vations were complemented with limited maps of the quiescent- 
gas tracers N2H + (l-0) and C I8 0(2-1). All observations were 
made with the dual-band, dual-polarization, single-pixel EMIR 
receiver. For the simultaneous N2H + and C 18 observations, we 
used frequency-switching mode to ensure that the spectra were 
free from ambient-cloud contamination. For the rest of the ob- 
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Fig. 1. Structure of the HH 1 14 MMS outflow from CO(2-l) data. Left: maps of intensity integrated in velocity ranges symmetrically separated 
from the ambient cloud velocity (V = -0.5 km s -1 from N 2 H + data). First contour and contour interval are 1 K km s -1 in the high-velocity map 
and 5 K km s -1 in the low-velocity map. The star symbol at the offset origin represents the position of HH 114 MMS (or(J2000) = 05 18 m 15f2, 
6(J2000) = +7°12'2"), the two squares are the positions where the line survey was carried out, and the dots are observed positions. The dashed 
green line is the direction of the position-velocity diagram shown in the next panel. Right: position-velocity diagram along the outflow axis. To 
emphasize the low-level outflow emission, the contours are spaced by factors of square root of 2 starting at 0.1 K. 



servations, we used the same wobbler-switching mode used in 
the CO(2-l) observations. As spectrometers, we used the FTS 
and VESPA instruments, which, depending on the line, provided 
velocity resolutions between 0.06 and 0.7 km s _1 . The resulting 
spectra were processed with CLASS as done with the CO(2-l) 
spectra. 

3. Physical properties of the HH 114 MMS outflow 

Fig. Q] summarizes the main features of the CO(2-l) data. The 
two maps in the left and middle panels show that the CO 
emission has a markedly bipolar distribution with respect to 
HH 1 14 MMS, with the redshifted emission toward the northeast 
and the blueshifted emission towa rd the southwest . This distri- 
bution agrees with that found by lArce & Sargentl (|2006) from 
interferometric observations of the central arcminute around the 
YSO, and with the distribution of blue CO emission mapped by 
iLee et al. (2002), also with an interferometer (and restricted to 
the southern blue lobe). Our more extended CO(2-l ) maps show 
that the HH 1 14 MMS outflow has a high degree of collimation 
at all velocities, and that its emission extends for at least 200" 
(0.4 pc for the assumed distance of 450 pc). 

To illustrate the velocity field of the gas, we present a 
position-velocity diagram along the outflow axis in the rightmost 
panel of Fig. Q] This diagram shows that both outflow lobes fol- 
low an approximate linear increase of the terminal velocity with 
distance from the central source. The "Hubble-law" velocity pat- 
tern, however, is only approximate, since the fastest velocities 
are not reached at the end of the lobes but at an intermediate dis- 
tance. This can also be seen in the velocity maps of Fig.Q] which 
show that the gas faster than 20 km s _I forms two bright spots 
about 50" from the YSO. The origin of this velocity pattern is 
unclear, and the 1-beam sampling of our data is too coarse to al- 
low further investigation. The blue-red symmetry of the pattern, 
however, suggests that it is intrinsic to the outflow, and not the 
result of an interaction with the ambient cloud. While unusual, 
the pattern seems not unique to HH 1 14 MMS. Similar cases, to- 
gether with their difficult interpretati on in terms of existing out - 
flow models, have been discussed bv lChernin & Massonl(ll995l) . 

We complete our analysis estimating the energetics of the 
moving gas traced by CO(2-l). We follow standard practice and 
do not correct for the uncertain inclination angle of the outflow. 
The lack of overlap between the blue and red lobes, together with 
the lack of color mixing inside each lobe, make HH 1 14 MMS a 



Case 2 outflow in the classification o flCabrit & Bertoutl (Il986l) . 
suggesting an inclination angle between 15° and 75° (we esti- 
mate an outflow opening angle max = 15° from the CO maps). 
As a result, the inclination correction to the velocity lies be- 
tween 4% and a factor of 4, and this makes our energetics esti- 
mate a strict lower limit. To avoid contamination from ambient- 
cloud material, we also limit our analysis to emission faster than 
3 km s _1 . We also assume that the CO(2-l) emission is optically 
thin and in local thermodynamic equilibrium (LTE), with an ex- 
citation temperature of 25 K (from the CO(3-2)/ CO(2-l) ratio. 
Sect. j 5]l, and use a CO abundance of 8.5 x 10" 5 dFrerking et al.l 
119821) . 

With these assumptions, we estimate a total outflow mass of 
0.15 M G , a momentum of 1.3 M km s _1 , and a kinetic energy 
of 3.5 x 10 44 erg. These values are between four times less than 
(mass and momentum) and equal to (energy) those estimated for 
L 1 157 by Bachille r et al.1 (1200 ll) . who do not correct for incli- 
nation in this also Case 2 outflow. 



4. Maps of dense-gas tracers 

While the CO maps of the HH 114 MMS outflow are not 
qualitatively different from those of oth er young outflows (e.g., 
ILee et al.l 120021 : lArce & Sargentl |2006), the maps of the high- 
density tracers are truly remarkable. Fig. |2]shows that the emis- 
sion of dense gas tracers such as HCN, CS, SO, CH3OH, or 
SiO is not concentrated toward the central source, but presents 
two bright peaks toward the outflow lobes. HCO + , while cen- 
trally concentrated, also shows an elongation in the direction of 
the outflow, and only N2H + , and C ls O (the last one not strictly 
a high-density tracer) present compact emission peaks towards 
HH 114 MMS. (The southern N 2 H + and C lg O peak coincides 
with the vicinity of IRAS 05155+0707.) 

The similarity between the emission of the dense-gas tracers 
and the CO outflow is not limited to the integrated maps of Fig. [2] 
Velocity maps using the intervals of Fig. Q] show the same blue- 
red bipolar distribution and the same offset between high- and 
low-velocity peaks. Individual spectra, like those in Fig. [3] show 
prominent wings of clear outflow origin. 

Bipolar outflows are known to sometimes enhance the emis- 
sion of certain molecular species, especially SiO and CH3OH, 
due to the disruption of dust grains and the ices that coat them 
(Ivan Dishoeck & Blake 1998). Few outflows, however, display 
the dramatic impact seen in the emission of other tracers associ- 
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Fig. 2. Emission from dense-gas tracers in the vicinity of HH 114 MMS. The maps represent emission integrated over the complete velocity 
range of detection. Contours are at 10, 20, ... 90% of the map maximum, whose value in K km s -1 is given within brackets under the name of the 
transition. Central positions, dots, and solid star as in Fig.[Tj The open star in the last two panels represents IRAS 05155+0707. 



ated with dense core gas, such as CS and HCN. These species are 
routi nely used to trace the innermost environme nt around YSOs 
(e.g JAfonso et al.lll998tlLaunhardt et all f998). and are not ex- 
pected to represent outflow gas. As Fig. |2]shows, this is not the 
case in HH 1 14 MMS, where the core emission is dwarfed by 
the outflow contribution. 

As mentioned before, only a small group of outflows present 
the prominent emission from dense-gas tracers found towards 
HH 114 MMS. The best-known member of this group is the 
L 1157 outflow, which has long been recognized for present- 
ing unusually strong emiss i on fr om shock-tracer molecules, 
such as SiO dMikami et all 119921) . The molecular maps of 
iBachiller et al.l (120011) (their Fig. 6) show that most dense-gas 
tracers in L 1157 (including HCN, CS, SO, CH 3 OH, and SiO) 
present a remarkable pattern of outflow-dominated emission, 
similar to that shown in Fig.|2]for HH 1 14 MMS. 

The dra matic effect of the L 1 157 outflow in the chemistry 
of its gas led Bachiller et al. (2001) to coin the term "chemically 
active" to refer to t his outflow and others with similar prop erties, 
such as BHR 71 ( iBourke et aLlll997h iGarav et al.lll998l) . The 
maps of Fig. |2]show now that the term can also be applied to the 
HH 1 14 MMS outflow, which appears as one of the finest ex- 
amples of this class, thanks to its clean pattern of a single bright 
emission peak near the end of two lobes of similar size. 



5. Column density and abundance estimates 

To further characterize the chemistry of the HH 1 14 MMS out- 
flow, we convert the observed intensities into molecular column 
densities. For this, we use the results of a small survey of transi- 
tions at 1, 2, and 3 mm wavelength toward (60", 40") and (-50", 
-30"). Representative spectra from this survey are shown in 
Fig. H The full list of detections include CO(2-l), CO(3-2), 
CS(2-1), CS(3-2), SiO(2-l), SiO(3-2), SO(3 2 -2i), SO(4 3 -3 2 ), 
CH 3 OH(2 jt -l jt ), CH 3 OH(3 i -2 jt ), HCN(l-O), HCO + (1-0), and 
H 2 CO(2 02 -l i). 



At the time of the observations, and based on our CO maps, 
the two survey positions were thought to represent the bright- 
est peaks of each outflow lobe. Later mapping in multiple 
molecules, however, showed that these positions were offset 
from the true emission peaks by ~ 15". While these small off- 
sets are not critical for the abundance results, which seem al- 
most position-independent, they may affect the comparison of 
the HH 1 14 MMS outflow intensities with those of other objects. 

As a first step in our column density estimate we chose for 
each survey position a 20 km s _1 -wide interval of outflow veloci- 
ties (see values in Table[T]l. The molecular emission in this inter- 
v al was analyzed w ith an LTE method similar to that presented 
in lTafalla et al.1 (l2010h . although no dilution factors were applied 
to the data because the observations were not centered on the 
emission peak. For the CO emission, we used data from 5-point 
crosses with a 5" sampling to convolve the CO(3-2) emission 
to the 12" resolution of the CO(2-l) data. These two intensities 
were used to derive CO excitation temperatures between 25 and 
30 K and the column density values reported in Table Q] For the 
dense-gas tracers with two observed transitions (CS, SiO, SO, 
and CH 3 OH), we used population diagrams to determine simul- 
taneously column densities and excitation temperatures, which 
range from about 4 to 8 K and indicate that the molecules are 
subthermally excited. Finally, for the three species with only 
one observed transition (HCN, HCO + , and para-H 2 CO), we as- 
sumed an excitation temperature of 6 K (mean value of CS, SiO, 
SO, and CH 3 OH), and used the integrated intensities to estimate 
molecular column densities. 

Using the previous column densities, we determined for 
each species an abundance with respect to CO, and compared 
this value with the typica l undep leted value in a dense core as 
discussed in Tafall a et al. I J20Tol) . The results of this estimate 
are presented in Table [T] in terms of the enhancement factor 
/enh> which is the ratio between the outflow and the core CO- 
normalized abundances. As can be seen, the / en h factors in the 
two outflow positions agree with each other within a factor of 2 
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Fig. 3. Sample spectra from the line survey carried out toward two 
outflow positions. The offset of each position is indicated in the top 
panel (in arcseconds and referred to the map center of Fig. [TJ. The 
vertical dashed line indicates the velocity of the ambient cloud. 



for all observed species. This agreement, which is probably at 
the level of our measurement uncertainty, suggests that the two 
outflow lobes are chemically very similar. 

The results in Table[T]also show that the enhancement factors 
of the different species span a wide range of values. The small- 
est factor is that of HCO + , and is consistent with this molecule 
having little or no enhancement in the outflow gas. This result 
is consistent with the map in Fig. [2] which shows no emission 
peak for this molecule toward neither of the outflow lobes (al- 
though the emission is elongated in the outflow direction). All 
other species in Table Q] present significant abundance enhance- 
ments, and range from about one order of magnitude in the case 
of HCN and CS to more than three orders of magnitude in the 
caseofCH 3 OHandSiO. 

While large, the enhancement factors in Table Q] are com- 
parable to those in other young outflows. This can be seen 
by comparing the values in the table with those presented by 
Tafal la et al.l d2010l) . who measured enhancement factors in the 
L 1448 and IRAS 04166 outflows and re-evaluated the L 1 157- 
Bl es timates previously made by Bachill er & Perez Gutierrez! 
d!997l) . The values for these three outflows are in good agree- 
ment with our HH 1 14 MMS estimates, and for example, the 
geometrical mean of the / en h estimates in the two HH 1 14 MMS 
survey positions agrees with those in L 1157-B1 better than a 
factor of 3 for all molecular species. Such an agreement in en- 



Table 1. Survey column densities and abundance enhancements. 



(60", 40") a (-50", -30") b 



MOL 


Wfrtn _2 i 

1 V ^111 J 


J enh 


1 V ^111 J 


J enh 


CO 


1.8 10 16 




2.1 10 16 




HCO + 


1.0 10 12 


1.1 


1.2 10 12 


1.2 


HCN 


2.1 10 13 


11 


2.1 10 13 


7.2 


CS 


1.8 10 13 


14 


2.8 10 13 


19 


p-H 2 CO 


8.2 10 12 


68 


9.5 10 12 


65 


SO 


3.4 10 13 


290 


4.4 10 13 


320 


CH 3 OH 


4.6 10 14 


1.8 10 3 


3.2 10 14 


1.1 10 3 


SiO 


5.5 10 12 


5.3 10 3 


1.4 10 13 


1.1 10 4 



Notes. <a) Quantities measured in the Vlsr range 4 to 24 km s . 
(b) Quantities measured in the Vlsr range -26 to -6 km s -1 . (c) Abun dance 
enhancement factor over dense core value as in Tafalla et ahl feOlOh . 

hancement factors confirms the idea that HH 114 MMS must be 
a bona fide chemically active outflow. 

Although the abundances with respect to CO in the 
HH 114 MMS outflow are comparable to those in L 1157, the 
integrated intensities towards its brightest molecular peak are 
lower than those of L 1 157-B1 by an average factor of 2.5. This 
suggests that the HH 1 14 MMS outflow may not have accumu- 
lated as much enriched material as L 1157, although beam di- 
lution may still affect the HH 114 MMS intensities. The more 
southern location of HH 1 14 MMS, on the other hand, makes this 
outflow a good candidate for ALMA observations, something 
impossible with L 1 157 due to its high declination. The newly- 
identified chemically active HH 1 14 MMS outflow therefore has 
significant potential to shed new light on the still-mysterious na- 
ture of this outflow class. 
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